Algorithms for rank-de cient and ill-conditioned Toeplitz least-squares and QR factorization

نویسندگان

  • S. Thirumalai
  • K. A. Gallivan
  • P. Van Dooren
چکیده

In this paper we present two algorithms-one to compute the QR factorization of nearly rank-deecient Toeplitz and block Toeplitz matrices and the other to compute the solution of a severely ill-conditioned Toeplitz least-squares problem. The rst algorithm is based on adapting the generalized Schur algorithm to Cauchy-like matrices and has some rank-revealing capability. The other algorithm is based on adapting the augmented systems method to Toeplitz matrices. This algorithm does not suuer from the numerical inaccuracy of most Levinson-and Schur-based schemes proposed in the literature thus far because it avoids forming the normal equations either implicitly or explicitly. Some comparisons with more recent Toeplitz least-squares algorithms are also presented.

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Parallel Algorithms for Toeplitz Systems

We describe some parallel algorithms for the solution of Toeplitz linear systems and Toeplitz least squares problems. First we consider the parallel implementation of the Bareiss algorithm (which is based on the classical Schur algorithm). The alternative Levinson algorithm is less suited to parallel implementation because it involves inner products. The Bareiss algorithm computes the LU factor...

متن کامل

A weakly stable algorithm for general Toeplitz systems

We show that a fast algorithm for the QR factorization of a Toeplitz or Hankel matrix A is weakly stable in the sense that RR is close to AA. Thus, when the algorithm is used to solve the semi-normal equations RRx = A b, we obtain a weakly stable method for the solution of a nonsingular Toeplitz or Hankel linear systemAx = b. The algorithm also applies to the solution of the full-rank Toeplitz ...

متن کامل

Superfast and Stable Structured Solvers for Toeplitz Least Squares via Randomized Sampling

We present some superfast (O((m + n) log(m + n)) complexity) and stable structured direct solvers for m× n Toeplitz least squares problems. Based on the displacement equation, a Toeplitz matrix T is first transformed into a Cauchy-like matrix C, which can be shown to have small off-diagonal numerical ranks when the diagonal blocks are rectangular. We generalize standard hierarchically semisepar...

متن کامل

New Fast Algorithms for Structured Linear Least Squares Problems

We present new fast algorithms for solving the Toeplitz and the Toeplitz-plus-Hankel least squares problems. These algorithms are based on a new fast algorithm for solving the Cauchy-like least squares problem. We perform an error analysis and provide conditions under which these algorithms are numerically stable. We also develop implementation techniques that signiicantly reduce the execution ...

متن کامل

Fast Qr Factorization of Low-rank Changes of Vandermonde-like Matrices

VANDERMONDE-LIKE MATRICES LUCA GEMIGNANI Abstract. This paper is concerned with the solution of linear systems with coe cient matrices which are Vandermonde-like matrices modi ed by adding low-rank corrections. Hereafter we refer to these matrices as to modi ed Vandermonde-like matrices. The solution of modi ed Vandermondelike linear systems arises in the approximation theory both when we use R...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

عنوان ژورنال:

دوره   شماره 

صفحات  -

تاریخ انتشار 2007